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H ' Abstract 

'^ , We present high statistics data on the distribution of shortest path lengths between two 

^ I near-by points on the same cluster at the percolation threshold. Our data are based on 

Q ' a new and very efficient algorithm. For d = 2 they clearly disprove a recent conjecture 

by M. Porto et al, Phys. Rev. E 58, R5205 (1998). Our data also provide upper 
bounds on the probability that two near-by points are on different infinite clusters. 
> 

cn : 1 Introduction 

^ . Although percolation is a problem which has been studied in great detail during the last 

G^ ! few decades fll, and although many exact results are known by now, there are still open 

questions which either have not been studied at all, or which have not yet been understood. 

In the present note we study spreading of percolation in form of an epidemic process 

0. In the physics literature, this is often called Leath growth of clusters @]. In this process 

^ . (which we assume to proceed in discrete time steps) one starts with an infected seed as single 

O '. 'growth site', and keeps at each time a list of growth sites. In the next time step, the list 

^ ' consists of all wettable sites which are nearest neighbors to any of the present growth sites, 

while the old list is canceled. A site is wettable if it had not been a growth site before, and 

^ ■ if it can be occupied (in site percolation) or it is connected to the present growth site by an 

occupied bond (in bond percolation). It is well known that the distribution of growth sites 

satisfies at the critical point and for large times t the scaling law P, Q 



Here, /3 is the well known critical exponent governing the fraction of sites occupied by the 
infinite cluster [jl| , z/^ is the critical exponent governing the correlation time if one goes off- 
critical, and z = f /ft is a dynamical exponent. The inverse of the latter is often called dmin- 
Finally, 0(C) is a universal scaling function. In a purely geometrical interpretation of eq.(l), 
t is often called the 'chemical distance' between the site x and the origin, as it counts the 



number of lattice steps on the cluster needed to reach x from the origin (i.e. the length of 
the shortest 'chemical path'). 

The main problem we want to study is the behavior of the scaling function (f){() for small 
(. This has been studied recently by Porto et al. (see also 0). As usual one expects a 
power law, 

m-c'. c-o. (2) 

The authors of used an analogy with self avoiding random walks (SAW) to conjecture 

Qi = c^min — /^/^^ (conjectured) (3) 

where d is the dimension of the lattice in which the cluster is embedded. This conjecture 
was then checked numerically and found to be satisfied [] 

There are however a number of problems associated with that analysis. The first is that 
the analogy with SAW is not very stringent. The analogous equation to eq.(|]) had been 
derived by deGennes using the fact that the end point density is the order parameter field 
in the SAW problem. Therefore, there is no anomalous critical exponent for the number of 
self avoiding loops (which have no end points), and this then gives an expression for the 
probability that a SAW nearly forms a loop which formally seems to involve mean field 
exponents. The authors of [§ essentially used these 'mean field' arguments, although there 
is no analogous underlying field theory (the field theory for percolation is Potts-like |]^), and 
we see no logical basis for eq.(0). 

Similarly, the simulations used in |]^ seem far from conclusive. In these simulations, 
clusters have been studied with t < 1800 in d = 2 and with t < 800 in d = 3. In both cases, 
'more than 100,000' clusters have been analyzed, at several fixed values of t. In figure 1 we 
show results for d = 2 obtained by their method. Each curve is based on more than 3 x 10^ 
clusters. In addition, in order to enhance the statistics, we have lumped together data with 
0.9t < t' < t in each curve. In this way we arrive at statistics at least 10^ larger than those 
of 0. In order to see more details we do not show simple scaling laws as in [^, but show 
data multiplied by a suitable power of r/t^. If eq.(H) holds, we expect to see a horizontal 
line for small x-values. Although our data are certainly not in contradiction with this, the 
very large statistical and systematic deviations from such a line render any precise statement 
impossible. 

If we want a significant test of eq.(|), we have thus to proceed differently. Indeed, for 
fixed r there is a much faster numerical method. While the above simulations need a CPU 
time ex t^'® in order to analyze one cluster, our improved method is faster by roughly one 
power of t. 

Let us assume we want to test whether two sites x and y on a regular (hypercubic) 
lattice are on the same cluster, and want to measure the length of the shortest connecting 
path if they are. Let us assume furthermore that the distance between x and y, measured 



""^In the present paper, p(x, t) is the density averaged over all clusters containing the origin, not only over 
the infinite one. It seems that the authors of ||5|| wanted their argument to apply to the incipient infinite 
cluster only. They implemented this by analysing only clusters with chemical radius tmax > 2000, from 
which they presented results for chemical distances up to 1800. If there were a substantial dependence on 
the chemical radius, this would lead to spurious violations of scaling since t,nax introduces a new scale. But 
this dependence seems to be weak, as also assumed in H. 
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Figure 1: Radial density of sites infected at time t, with arbitrary normalization, plotted against 
r/t^ with z = 0.8844 (corresponding to dmin = 1.1307 it 0.0004; ^M). In order to enhance the 
significance of the plot, each cm've corresponds to times e]0.9t,t], and the data are divided by 
rrfmi„-/3/v. The curves correspond to t = 60, 100, 140, 200, 300, 500, 800, 1200 and 1800. 



as the sum of the coordinate distances, is even (the case of odd distances will be discussed 
below). Instead of growing a single cluster from a seed at either x or y, we grow two clusters 
simultaneously, one starting from x and the other from y. We stop the growth in any of the 
following cases: 

(1) Both clusters have a common growth site. If this happen at time t, then there exists 
a path of length 2t passing through this growth site and connecting x with y. Since this is 
the shortest such path (otherwise the growth would have stopped before), in this case the 
chemical distance between x and y is 2t. 

(2) The cluster growing from x dies (i.e., has no more growth sites). In that case, there 
cannot be a path from x which is long enough to reach y, and x and y are on different 
clusters. Notice that both clusters cannot overlap. They cannot ovelap in a site which has 
the same chemical distance from both seeds, because they would then have been killed by 
rule (1) above. And they cannot overlap in a site which, say, has a smaller chemical distance 
from X than from y. Such a site would first have been wetted by x. In order to be wetted also 
from y, it must have at least one other wettable neighbor which was not wetted, however, 
at the time step following its first wetting. This is not possible. 

(3) The cluster growing from y dies. 

(4) t reaches some upper bound T specified at the beginning. In this case x and y are 
either on different clusters, or the shortest path has length > 2T. 

If the distance between x and y is odd, we have to replace (1) by: 

(r) A grow site of the first cluster at time t coincides with a growth site of the other 

cluster at time t — 1. In this case, the chemical distance is 2t — 1. 

Let us denote by p{t) the probability that an event dies because of rule (1). It is equal 

to the probability that two sites at a distance r = |x — y| have chemical distance 2t, which 
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Figure 2: Histograms of n{t) and —dN{t)/dt for d = 2 bond percolation with x — y = (1, 1). In 
order to reduce statistical fluctuations, data are averaged over intervals [0.98t,f]. 



according to eqs.(l-3) is given by 

pit) 
with 



P(x-y,2t) 
E?LoP(x-y,t') 
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Inserting here the conjectured eq.(|3), we arrive at 



A = 2 + 



^ 
vt 



(conjectured) 



(4) 
(5) 

(6) 



During these simulations, we not only collect a histogram n(t) of the times t when rule 
(1) applies (which gives the distribution p(t) when normalized), but also a histogram of times 
when the growth stops due to any reason or, equivalently, a histogram A^(t) indicating how 
many events have survived at least t time steps. The numer N{T) of events surviving until 
the very end gives an upper bound on the probability that x and y are on different infinite 
clusters, Pdis < N(T)/N. If Pdifr = (which we expect for d < 6; there are several infinite 
clusters for any d > 2, but the chance that they come close to each other should be zero), 
we expect N{t) also to decay with a power law. 



N{t) ~ t~ 



(7) 



with exponent ;U < A — 1. This should be compared to the power with which single clusters 
survive, —(3/vt [ii i|- Thus according to the conjecture the difference in the powers is exactly 
1, suggesting that the present algorithm is faster by one power of t. 

In fig. 2 we show, for d = 2 and x — y = (1, 1), both n{t) and \dN{t)/dt\. We used bond 
percolation where Pc = 1/2 exactly. We see large deviations from power laws, but these 
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Figure 3: Plot of N{t)t^'^^^^ for d = 2 with x— y = (1, 0) (lower curve) and for x — y = (1, 1) (upper 
curve). The numbers of configurations were chosen such that both curves coincide for large t. The 
exponent 1.0921 = 1 + f3/ut is such that the curves should become flat for large t, according to the 
conjecture of [^. Relative statistical errors are ~ l/^yN{t). For t = 4000, they are ~ 2.6 x 10~^. 
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Figure 4: Similar to fig.3, but for d = 3. This time we used f3/ut = 0.3467 ± 0.0028 from Q, 
Relative errors for N{t = 800) are ^ 3 x 10~^. 



deviations are the same for both curves. Obviously, if a pair of clusters dies, the chance that 
it dies because of rule (1) is finite and tends to a constant, and therefore 

/i = A - 1 . (8) 

If we accept this, we can obtain the most precise estimate of A from N{t). In any case, 
even if this is not correct, we obtain from N{t) a lower bound on A. Since we shall see that 
this lower bound is larger than the value conjectured by 0, this is sufficient to exclude the 
conjecture. 

Results for N{t) for rf = 2 are shown in fig. 3, for x — y = (1,1) and for x — y = (1, 0). 
Each curve is based on > 10^ runs with T = 4000, and took about 90 h CPU time on 

a fast workstation. In order to compare with eq.(^), we multiplied N{t) by t '^, using 
the exact value /? = 5/36 [l[| and the estimate z/j = 1.5075 ± 0.0004 (corresponding to 
c^min = 1.1306 ± 0.0003). This estimate is based on new simulations which follow exactly 
the lines of [0, but have four times higher statistics. It is fully compatible with [|I^. We 
see clearly that eq.(P) is wrong since it gives a too small value. A precise estimate of A is 
hampered by the very strong corrections to scaling visible in fig. 3. Fitting these corrections 
by terms ~ 1/t relative to the leading term and assuming that n(t) oc N{t), we arrive at 

/i = 1.1055 ±0.0010, A = 2.1055 ±0.0010, ^i = 1.041 ± 0.001. (9) 

Basing the analysis on n{t) instead of A^(t), we would get A = 2.105±0.002, gi = 1.041±0.003. 
A similar analysis for d = 3 is shown in fig. 4. This time T = 800. For pc we used the 
value 0.2488126 of 0. The most precise value of P/ut can be obtained either by combining 
the value of z from |^ with the value of r from |^, or by using directly the estimate of (3/i't 
from §. The first gives (3/ut = 0.347 ±0.004, the latter p/ut = 0.345 ±0.004. In addition we 
performed further simulations using the method of 0. Together, all these combine to our 
final estimate jS/ut = 0.3467 ±0.0028. As seen from fig. 4, the corrections to scaling are even 
worse than for d = 2, and the disagreement with eq.(^ is much less pronounced. It seems 
that eq.(P) is just barely compatible with our data which give, after taking into account the 
corrections to scaling, 

/i = 1.353 ± 0.003 , ^1 = 0.905 ±0.008. (10) 

We also tried to fit N{t) by an exponential plus a constant, the latter corresponding to 
events where the two starting points are on different infinite clusters. We obtained Pjiff < 
10"^ for rf = 2 and Pdiff < 2 x 10"^ for d = 3, both for x - y = (1, 0) and for x - y = (1, 1). 

Finally, we should point out that the advantage of our two-seed algorithm over the naive 
one using a single seed is even more pronounced for supercritical percolation. There, the 
survival chance for an event decays exponentially with t, while it would not decay at all for 
the growth of a single cluster. 
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